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1. Introduction 

Domain walls are two-dimensional membranes that form when a discrete symmetry is 
broken at a phase transition, e.g., the interfaces (Bloch, Ising walls) between magnetic 
domains in ferromagnetic materials. In Cosmological models of the early Universe domain 
walls form according to a pattern known as the (second) Kibble mechanism; see [1]. The 
impact of domain wall defects and of other topological defects in Cosmology has been 
studied in depth in the monograph [2] . The evolution of domain wall networks is a problem 
of particular interest in this context, see e.g. [3]. This problem is usually studied in 
computer simulations, although an analytic approach has been developed in [4] and [5]. 
By identifying the moduli space of domain wall networks in a U{Nq) gauge theory with 
Nf scalar fields in the adjoint representation, the authors of [4] and [5] implemented the 
low-energy dynamics of the network by studying the geodesic motion in the domain wall 
network moduli space. It is also interesting to consider domain walls as the seed of Randall- 
Sundrum scenarios [6], in which space-time is five-dimensional and the 3-brane is wrapped 
around some background four-dimensional gravitational field, while the particle dynamics 
is concentrated inside the wall; see [7]-[8]-[9]-[10]-[ll]. In this framework some authors have 
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considered the possibility that the Big Bang of the standard (3+1) dimensional cosmology 
was originated from the collision of two branes within a higher dimensional spacetime, 
leading to the production of a large amount of entropy [12]. 

Our purpose here is to investigate the very rich moduli space of domain walls in a 
hybrid of the non-linear sigma model and the Ginzburg-Landau theory of phase transitions. 
The linear Q(iV)-sigma model is the key to our present understanding about the origin of 
mass: see [13] and [14]. The non-linear S^-model version only describes the dynamics of 
Goldstone bosons. In [15], we addressed an especially simple case of massive non-linear 
sigma model: we chose S 2 as the target manifold; assigned different masses to the two 
pseudo- Goldstone bosons, and were able to identify all the domain walls of the system 

Regarding [15], topological defects 1 in massive non-linear sigma models have been 
known for some time and have been profusely studied in different supersymmetric models 
under the circumstance that all masses of the pseudo-Nambu- Goldstone particles are equal. 
The study started with two papers by Abraham and Townsend [16], [17], in which the 
authors discovered a family of Q-kinks in a (l+l)-dimensional Af = (4,4) supersymmetric 
non-linear sigma model with a hyper-Kahler Gibbons-Hawking instanton as the target 
space and mass terms obtained from dimensional reduction. In [18], however, these kinks 
were re-considered by constructing the dimensionally reduced supersymmetric model by 
means of the mathematically elegant technique of hyper-Kahler quotients. By doing so, 
the authors dealt with massive CP N or MP N models: a playground closer to our simpler 
massive S 2 -sigma model. Similar Af — 2 BPS walls in the CP 1 -model with twisted mass 
were described in [19]. In a parallel development in the (2+l)-dimensional version of these 
models, two-dimensional Q-lumps were discovered in [20] and [21]. Within this field, the 
most interesting result is the demonstration in [22] and [23] that composite solitons in 
d = 3 + 1 of Q-strings and domain walls are exact BPS solutions that preserve \ of the 
supersymmetries: ( See also the review [24], where a summary of these supersymmetric 
topological solitons is offered.) 

Our research differed from the above works in two important aspects: 1) We stuck 
to a purely bosonic framework. 2) We studied the case when the masses of the pseudo 
Nambu-Goldstone bosons are different, a property that forbids extended supersymmetries. 
The search for domain walls in the d — 3 + 1-model is tantamount to the search for 
finite action trajectories in the repulsive Neumann system [25]: a particle moving in an 
S 2 -sphere under the action of non-isotropic repulsive elastic forces. It is well known that 
this dynamical system is completely integrable [26], [27]. We showed, however, that the 
problem is Hamilton- Jacobi separable by using elliptic coordinates in the sphere. Use of this 
property allowed us to find four families of homoclinic trajectories starting and ending at 
one of the poles which are unstable points of the mechanical system. In the field-theoretical 
model, the poles become ground states, whereas the homoclinic trajectories correspond to 
four families of non-topological domain walls. Each member in a family is formed by a 
non-linear combination of two basic topological domain walls (of different type), with their 
centers located at any relative distance with respect each other. 

1 Oi varying character, depending on the spatial dimension and the charge: kinks, Q-kinks, lumps, 
strings, walls, etcetera . 
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Here we shall address a Ginzburg-Landau non-linear S 2 -sigma model; i.e., we will keep 
the target space but we add a quartic, rather than quadratic, independent of field gradients, 
contribution to the potential energy density; see [28] for a mathematical definition of these 
models. Because of the constraint, the GL function must be non-isotropic and we shall 
consider non-equal quadratic and quartic couplings in such a way that the anisotropy is 
maximal. The consequence is the existence of a spontaneously broken discrete symmetry: 
a necessary condition for the existence of domain walls. 

We shall further restrict (but not too much) the space of parameters of the model, the 
quadratic and quartic couplings besides the radius of the § 2 -sphere, bearing in mind that 
we have to deal with an integrable analogous mechanical problem. Instead of the Neum- 
mann system we must solve the problem of a particle moving on the S 2 -sphere under the 
action of non-isotropic inelastic (non-harmonic) forces. Fortunately, this mechanical sys- 
tem is also Hamilton- Jacobi separable, and we shall apply the Hamilton- Jacobi procedure 
to find all the finite mechanical action trajectories by using sphero-conical coordinates. 
There are four unstable points -four ground states in the field theory-; three types of basic 
heteroclinic trajectories -topological domain walls- joining unstable points; four families of 
heteroclinic trajectories joining antipodal unstable points -topological domain walls- with 
the same mechanical action -wall tension-, and four families of homoclinic trajectories - 
non-topological domain walls- also with the same mechanical action 2 . We remark that the 
basic domain walls are usual walls, concentrated at a point. Domain walls belonging to 
any of the degenerate families are composite domain walls in the sense that, generically, 
the walls are centered at two points, resembling a non-linear superposition of two basic 
walls. 

In order to describe all this, we shall organize the paper as follows: In Section §. 
2 we introduce the model, explain the physical content, and describe the basic domain 
walls that can be found by applying the trial orbit method. Section §. 3 is devoted 
to addressing the analogous mechanical problem. Sphero-conical coordinates are used to 
show the Hamilton-Jacobi separability. A new basic wall is easily guessed, and the other 
two, previously known, are also expressed in these coordinates. In Section §. 4 we apply 
the Hamilton-Jacobi procedure in full generality. We find the families of non-topological 
and topological domain walls in an explicit analytic form. This is remarkable: the solution 
for the orbit and time-schedule equations provided by the HJ prescription is frequently 
expressed in an implicit form that is difficult to invert. We have succeeded, however, in 
performing the inversion in this problem. We offer a last Section §. 5 with further comments 
and some suggestions for future lines of enquiry. 

Finally we have complemented this paper with a MATHEMATICA file, which can be 
found at http://campus.usal. es/^mpg/General/Mathematicatools.htm. This file includes 
animated figures which display the behaviour of the domain wall families depending on the 
different coupling constants and family parameters. 



2 The wall tensions of the topological and non-topological families are different. 
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2. The Ginzburg-Landau non-linear 8 2 -sigma model 



The action and the constraint governing the dynamics of this hybrid of Ginzburg-Landau 
and non-linear Sigma models are respectively: 

Owing to the constraint, the fields take values in the § 2 -sphere of radius mR embedded in 
R 3 . Setting an ortho-normal frame, e ( a • e& = a, b = 1, 2, 3, in R 3 , we write the fields in 
the form: 

3 

£ : R i,3 _^ §2 . x(y») = Xa{y^)e a . 

a=l 

xiy^) are maps from the Minkowski space R 1 ' 3 into S 2 C R 3 . The contra- variant tetra- 
vector y^, \i = 0, 1, 2, 3, provides local coordinates for a point in R 1 ' 3 : the Minkowski space 
equipped with the metric tensor g^ v = diag(l, —1, —1, —1). Thus, y^y^ = y\ — y\ — y\ — y\, 

dy» dy» ~ dyl dy\ dy* dy\ ' etCetem - 

Throughout the paper we shall use the natural system of units, in which the Planck 
constant and the speed of light in vacuum are in the units: H = c = 1. Therefore, the 
physical dimensions of the fields and the parameters m and /3 a are those of inverse length, 
[Xa] — [fia] = [m] = L~ x , whereas a a and R are non-dimensional couplings. In terms 
of non-dimensional fields (p a = ^Xa, space-time coordinates x M = m/, and quadratic 
couplings ?7 2 = ^2/S 2 , the action and the constraint read 3 : 

s$ = Jd*x^d ll $-£r$-v(<t> 1 ,<h,<h)} , 4(^) + 4(^) + 4(^) = R 2 

(3 \ 2 3 

a=l J 6=1 

We shall address the maximally anisotropic model and, with no loss of generality, choose: 
a\ > ol\ > al > 0. 

The static homogeneous configurations for which the action is extremal are the critical 
points of V complying with the constraint: 

dV 

—- + 2A0 a = O , (2.1) 

where A is the Lagrange multiplier forcing the constraint. There is an important parameter 
in the system: 



3 In V ((/>!, 02, 03 ) we have dropped the irrelevant constant: R 2 ^ Vs 2 Vl + (1 — a 2 R 2 ){a\ — a§)^ 



2 
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If S 2 E (0, 1), equation (2.1) is solved by 18 critical points, but the following four points: 

4>l = R 2 5 2 , 0l = , 4>l = R 2 (l-5 2 ) = R 2 5 2 (2.2) 

are the absolute minima of V in § 2 . The rest of the critical points are maxima or saddle 
points. In fact, we will choose 6 in the < 5 < 1 range, because the number of minima 
is maximized and this circumstance provides a richer structure for the domain wall space. 
The Lagrange multiplier at the four minima (2.2) can be easily computed: 

A = <4vi - aivi 

2{a 2 - a 2 ) 

2.1 Solving the constraint: particle masses 

In order to show explicitly the physical content of the model -symmetry breaking pattern, 
particle masses, interaction terms (with and without derivatives) as well as the physical 
characteristics of the parameters- it is convenient to solve the constraint by choosing, e.g., 
01 and 02 as independent fields: §\ — R 2 — 2 — 2 . On 8 2 , the unconstrained action 
becomes: 

ff s2 [0i,0 2 ] = / o 1 2\<2 I ^ dt {^§ 2 (d/i01,<9 M 2 , 01, h) ~ ^§2(01,02)} 
D&idpfa 0) = - ^01^01 + ^02^02 + ~ ^_02_02 J 

%(0i,0 2 ) = i(0? + a 2 2 - «) 2 + ^0 2 + ^0 2 (2.3) 



Again, here we have dropped an irrelevant constant and redefined: — »• ^ 2 ^ 2 x M . The 
new parameters are defined in terms of the old ones: 



^3 



2 qj-qj ^-% 2 « = ^-^1 

^ 2 2 ' 2 2' ' / 2 2\9 ' ^ / 2 2\9 

The range of a is due to the inequalities between the a a 's. 

We also restrict the number of independent parameters, for reasons to be explained 
later, and set 4 : /3 = a 2 (j + a 2 R 2 ). In a last move, we reshuffle the forgotten additive 
constants in such a way that the function Vg2 takes the value of zero at its minima. The 
potential depends only on <r, 5 and R ( we shall denote always a 2 = 1 — a 2 ) 

V §2 (0i,0 2 ) = 1 (0 2 + a 2 ^ - ^2)2 + Ii?Va 2 ^ . (2.4) 
The set M of zeroes of V§2, see Figure 1, 

M = {vi = (i2(5, 0,i?5),^ 2 = (-125, 0,i?5),^ 3 = (-125, 0,-/25),v 4 = (RS, 0,-125)} , 

encompasses the four ground states of the model (centered at the constant classical solu- 
tions) at the four degenerate absolute minima of Vg2(0i,02). 



Alternatively: (a? - a|)(ry 2 2 - ryf) = R\a\ - a\) + - rg) . 
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Expanding V§2 around any of the vacuum points 
V&{R5 + <h.,<h) = 25 2 R 2 <t>1 + ^R 2 a 2 a 2 <t> 2 2 + 

+ 2srm4 + v 2 4) + \(<t>l + ° 2 4>l? 

we see that: (a) The particle masses are: n\ — 45 2 R 2 , 
ji 2 = cr 2 a 2 R 2 . (b) The Z2 x Z2 symmetry engendered 
by 0i — 0i and 02 — > — 02 is spontaneously broken 
to the last Z2 subgroup by the choice of vacuum, (c) 

There are two trivalent vertices with couplings 2SR 
Figure 1: Ground states plotted for and 2§R(J 2 (d) There are three tetravalent vertices 
S = and ii = 1. These values CQU j m g l 2 an( j cr 4 

for R and (5 will be maintained in all W1 cou P m S s 2 ' a ' an 2 

figures throughout the paper. The r61e of ^ as a coupling constant comes from 

the following expansion of the 9^03^03 term in the 

Lagrangian: 

' ^lTTT 7 + <P27^— 7 <Pl o ^027^ + 



^R2 ^2 _ ^2 _ ^2 #2 ^ ^ y ^ ^ 

i.e., an infinite number of vertices with two field derivatives arise proportional to powers 
of -^2 due to the geometry of the system. 

2.2 Solitonic domain walls 

The non-linear field equations of the system are: 

3 2 a - V 2 a = -2a 2 J a a Ub " 1 j " V 2 a <Pa + A a , a = 1, 2, 3. (2.5) 

We temporarily return to keep 03 explicit. A is again the Lagrange multiplier in the 
equation (2.5). For any solution of (2.5) it can be shown to be 5 : 



d - 

(2.6) 

where V is the gradient in the spatial subspace of Minkowski space. 

Our main goal in this paper is to investigate the domain wall solutions in this model. 
Domain walls are non-singular solutions of the field equations (2.5) such that their energy 
density has a space-time dependence of the form: £ (x°, x 1 , x 2 , x 3 ) = £(x 1 — vx°), where v 



5 Multiply (2.5) by <j) a and sum the three equations. Also use that: Y^ a =i ^ad^a = to perform a 
partial integration. 
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is some velocity vector in the x 1 direction, and their energy functional: 

E[4>] = lim — g / dx 1 ( \dv4> • <9 O + \di$ • <9i0 + V (fa, fa, fa)) 
L^oo af — a% J \2 2 J 

L 2 f 

= lim — k I dx 1 £(x°,x 1 ) , 

L^oo af — a% J 

is proportional to the area L 2 of a normalizing square in the X2 — £3 plane. Therefore, 
these solutions will be domain walls or solitonic (thick) 2-branes orthogonal to the a? 1 -axis. 
The Lorentz invariance of the model implies that it suffices to know the x°-independent 
solutions fax 1 ) in order to obtain the domain walls of the model: fax®, x 1 ) — fax 1 — vx°). 
For static and x 2 -, x 3 - independent configurations the PDE system (2.5) becomes the 
following system of three ordinary differential equations: 

= -2a 2 a <j> a (j2 a Ub ~ lj " Vl<t>a +\<t>a, a = 1,2,3 (2.7) 
and the tension of the wall (2-brane) reduces to: 

The ODE system (2.7) can be interpreted as the Newton equations of a mechanical 
system, which we shall refer to as the analogous mechanical system to our field theoretical 
problem. Thus, the x 1 coordinate in R 3 will be identified with r: the mechanical time. The 
field configurations fa(x 1 ) will give the paths in § 2 , X a (r). The V (fa, fa, fa) function of 
the field theory will be minus the mechanical potential V(X\,X*±, X%). Finally, the domain 
wall tension ft will be interpreted as the mechanical action functional. We shall always 
use, however, the field theoretical notation, although the interpretation should be clear. It 
should be stressed that the mechanical potential is minus the function V. 

The finite wall tension (finite mechanical action) requirement is fulfilled if and only if 
the asymptotic conditions hold: 

lim ^- = , lim (fi G A4 . (2.8) 

Thus, the space of finite wall tension configurations 

4 

C = {(f) e Maps(M x R 2 , § 2 )/Maps(R 2 , point) : tl[4>] < +00} = |J dj 

is the union of sixteen disconnected sectors: Cij, labeled by the element of M reached by 
each configuration at x 1 — > — 00 and x 1 — > 00. If i ^ j, the finite tension walls will be 
termed as topological walls, whereas non-topological walls will be the solutions belonging 
to the Cu sectors. 

2.3 Trial orbits and two basic walls 

Before searching for general domain wall solutions of the ODE system (2.7), we shall show 
two particular ones by sticking to Rajaraman's trial orbit method [29]. 
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In the meridian 02 — 0, which we choose 
as a trial orbit, the constraint becomes <j)\ + 
03 = R 2 . The polar angle in this maximal 
circle solves the constraint 



4>\{x) — R sin#(x) 



03 (x) = Rcos0(x) 



We have written x = x 1 for simplicity and will 
maintain this convention in the rest of the pa- 
per. The ODE system (2.7) reduces to the sin- 
gle second-order OD equation: 



d 2 6(x 
dx 2 




■i — e 



Figure 2: Graphics of V§i(0). 



R 2 sm29(x) (sin 2 9(x)-S 2 ) 



The mechanical potential on the orbit is: 



(2.9) 



U(6) = -V s i(0) 



R 4 



(sin 2 <9( 



x 



The mechanical energy I provides a first-integral for (2.9): 



T i (dey 1 i fde 



i? 2 



2\ 2 



(sin 2 0(x) - <5 2 ) 



The critical points of E7 are: (1) minima: 9q — 0, the North pole, and n — ir, the South 
pole and Ok = | and = antipodal points in the equator. (2) maxima: + q — arcsinJ, 
9-0 — — arcsinJ, 6 +n — arcsin<5 + ir, 9- w — — arcsinJ + ir, respectively, the v\, v%, v$, and 
V4 minima of the field theory. The mechanical energy for all these maxima is I = 0, 
which must therefore be the value of the integration constant of (2.9) required to obtain 
solutions with finite wall tension. Therefore, the topological wall solutions correspond to 
the quadratures of 

f- = ±R (sm 2 9 - S 2 ) 

dx v 7 

which produce two types of analytical outcomes: 

1.- Polar Meridian Domain Walls (PMW): We shall denote this kind of solutions as 
#f 2 MW (x), 9™ w (x), e^ w (x) and ^ 3 MW (x), where the indexes stand for the asymptoti- 
cally connected vacua via the domain wall. For example, the orbit of the solution #f 2 MW 
corresponds to the piece of the 02 = meridian joining the v 1 and v 2 vacua, which crosses 
the North Pole. 6^ MW is the anti-wall of the previous solution while 6^ 4 MW and 6>£ 3 MW are 
similar domain walls confined to the South Polar Region. We find that 



6™ w (x) = arcsin 



S sinh[RSS(x — xq)] 
fcoah 2 [RS5(x-x )] ~ S 2 



+o) 
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and therefore 9™(x) = 0™ w (-x), €^0*0 = 0™ w + tt G (0-^0+*), and 0f 3 MW (x) = 
0™ w (— These topological walls live in the topological sectors C12, C34, C21, and C43. 
All of them have the same tension, fi(PMW) = ft(0i 2 ) = ^(02i) = ^(034) = ^(043), where 

i? 3 



f](PMW) 



at 



t x -a 3 

We write all these topological walls in a unified way in the original field variables 

(-l) ei i?5sinh [RSS(x - x )] 



[66- (1 -25 2 ) arccos 5] 



^PMW 



iPMW 
9 3 



(x;ei) 
(x;e 3 ) 



cosh 2 - x )] - 6 2 

(-l) e3 i?5cosh [i?(W(x-a )] 
cosh 2 [R66(x — xq)] — 6 2 



</>™(x) = 



ei,e 3 = 0, 1 



(2.10) 



For 63 = the walls belong to C12 (ei = 1) and C21 (ei = 0). For 63 = 1 the walls belong 
to C34 (ei = 1) and C43 (ei = 0). The (^-component has the form of a kink and the 
03-component is bell shaped. 



am 
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Figure 3: Field profiles (2.10) for e\ = = 63, xq = (left). Orbits (red curves) in C21, C12 and 
C43, C34 (right). 

2.- Tropical Meridian Domain Walls (TMW): We shall now denote these solutions 
as 6™ w (x), 6™ w (x), 0j 3 MW (x) and 0j 2 MW (^)- F °r example, 0^ 4 MW (x) connects from the 
vacuum v 1 to the v 4 vacuum crossing the equator of the sphere. We find that 



qTMW 

'41 



(x) 



arccos 



5 smh[R55(x — xq)] 
/ cosh 2 [R6 6 (x-x )}-5 2 



C w We(i^ + o) 



and therefore 6™ w (x) = 9™ w {-x), 6™ w (x) = 6™ w (x)+7r e (0 +7r ,0_ o ) nnd6™ w (x) = 
^32 MW (— x). These topological walls belong to the C41, C32, C14, and C23 sectors and their 
tension is O(TMW) = Q(8™ w ) = O(0™ w ) = ft(0j 3 MW ) = O(0™ w ) where 

i? 3 

fi(TMW) = — j t M + C 1 ~ 2(52 ) arccos 5 ] • 

^2 

In the original field coordinates the analytical expressions are: 



j TMW/ x (~1) K1 RS cosh [R88(x - x ) 



cosh 2 [R58(x - x )] - £ 2 



b™ w (x) = 
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(-l) K3 Msinh [RSS(x - x )] 
cosh 2 [R55(x — xq)] — 8 2 



Kl,K 3 = 0,1 



(2.11) 



If ^3 = 0, the domain walls live in C41 (fti = 0) and C14 (fti = 1). If ^3 = 1 the domain 
walls live in C32 (fti = 0) and C23 (fti = 1). The component is bell shaped and the ^3 
now has the form of a kink. 





Figure 4: Field profiles (2.11) for ft 1 = = ft 3 , xq = (left). Orbits (red curves) in C41, C14 and 
C23, C 32 (right). 

To end this subsection we show ( Figure 5 ) the tension densities of these two kinds of 
topological wall. 




Figure 5: Cross-sections of energy densities: (left) PMW walls in C12 and C21. (right) TMW walls 
in C23 and C14. 



3. Hamilton-Jacobi separability 
3.1 Sphero-conical coordinates 

In order to search for all the domain wall solutions of the model, we introduce sphero-conical 
coordinates (Ao, Ai, A2) in M 3 , see, e.g., [25]-[26]: 

,2 \ ( a i ~ Ai)(ai - A 2 ) 2 \ ( a 2 - Ai)(a 2 - A 2 ) 2 a ( a 3 ~ Ai)(a 3 - A 2 ) 

01 — A 07 w 7? 02 — A 7 w ^ 03 — A 07 w ^ 

(ai - a 2 )(ai - a 3 J (a 2 - ai)(a 2 - a 3 J (a 3 - ai)(a 3 - a 2 J 

with separation constants related to the cr-parameter: 

a>i = , a 2 = a 2 , a 3 = 1 4=> < Ai < a 2 < A 2 < 1 
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In this system of coordinates, the constraint is simply Ao = <j)\ + (f) 2 + (\)\ = P 2 , such that 
the field components restricted to S> 2 in terms of sphero-conical coordinates read: 



6? = — 
1 a 2 



Ai A2 



R 2 
a 2 a 2 



(a 2 - XtXXi - a 2 



R 2 



(1-Ai)(l-A 2 ) 



(3.1) 



The map induced by the change of coordinates (3.1) is eight-to-one; i.e., each octant of the 
S> 2 sphere is mapped onto the rectangle P2 in the (Ai, A2)-plane. See Figure 6 for a — -^=, 
a selection of a maintained in all the graphics below. 




1 
1 
1 
1 



Figure 6: The S 2 sphere in R 3 (left). The P 2 rectangle in the (Ai, A 2 )-plane (right). The dotted 
green line of P 2 is mapped to the equator while the dashed blue line of P 2 is mapped to the meridian 
which crosses the foci. The red solid line in P 2 corresponds to the other meridian displayed in the 
sphere. The Ai = constant and A 2 = constant iso-curves are shown back in S 2 . 



The sphero-conical coordinates distinguish four special points in 8 2 : the foci Pi, P2, 
P3 and P4, all of them mapped onto the corner (Ai,A2) = (a 2 , a 2 ) of P2. There is a 
direct relation between these coordinates and the elliptical coordinates on a sphere used in 
[15]- [30]. Choosing two non-antipodal foci, for instance Pi and P2, we have that: 



Ai — sin 



r\ - r 2 
2R 



Xo — sin 



2 ( n + r 2 
2R 



where r\ and r 2 are the geodesic distances from 
a given point in § 2 to Pi and P2 respectively. 



T\—T2 



and u 



2 



are the spherical-elliptic 



coordinates in S 2 . Thus the iso-curves depicted 
in Figure 6 (left) represent "ellipses" and "hyper- 
bolas" on § 2 . 

The eight-to-one correspondence between § 2 
and P2 maps all the Cij sectors onto only one; 
also, the four points ^1,^2,^3 and V4 are mapped 
onto the point v = (a 2 , 5 2 ) in P2. See Figure 7. 

We are assuming that a 2 < 5 2 . The choice 
a 2 > 5 2 is equivalent to this modulo a | rotation 
around the (/>2-axis. 



1 

h 

1 

-If 



Figure 7 : The ground states and the foci 
in sphero-conical coordinates. 
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It is easy to check that the non-derivative part of the field theoretical energy density, 
the ^-function (2.3), in sphero-conical coordinates reads: 

%(Ai, A 2 ) = — ((a 2 - X^S 2 - X 1 ) 2 + (A 2 - a 2 )(A 2 - 8 2 ) 2 ) 

if and only if j3 = a 2 (7 + a 2 R 2 ). 

The action in sphero-conical coordinates is 

S = f d 4 xC = J d 4 x |^ii(Ai,A2)^Ai^Ai + ^ 2 2(A 1 ,A 2 )a M A 2 ^A2-y s2 (A 1 ,A 2 )| 

where the components of the metric tensor induced in P2 by the change of coordinates are: 

ii a x , -1 4A 1 (a 2 -A 1 )(l-A 1 ) 22 , 4A 2 (A 2 -a 2 )(l-A 2 ) 
5 (^2)=^= ^ 2(A2 _ Ai) , 9 (Ai,A 2 )= 522 = ^ 2(A2 _ Ai) • 

Defining 

< = tI^-=5h(Ai,A 2 )^A 1 , ^ = ^= 522 (A l5 A 2 )^A 2 , 

one sees that the energy-momentum tensor = -K^d v \\ + 7r 2 9 !/ A 2 — g^/C of the wall 
solutions of (2.7) is diagonal, in the form: 

/ 3Ai5Ai <9A 2 <9A 2 \ nnin 



Tl(x) — is due to the continuity equation d x T^(x) = and the remaining components 

r;(x) oc s» 



are of the form T"(x) oc ^ because of the parallel unbroken Lorentz invariance (in the x 2 -, 



x 3 -directions) of the wall. 
3.2 The analogous mechanical system 



The mechanical momenta are: tt\ = gn(Ai, A2)^r, ^2 = Sn(Ai, ^2)^ • The mechanical 



Hamiltonian 

2A 1 (cx 2 -A 1 )(l-A 1 ) 2 2A 2 (A 2 -cx 2 )(l-A 2 ) 2 
H = i?(A 2 -A!) " 1 + i?(A 2 -A!) ^ + C/(Ai,A 2 ) , 

f^(Ai, A2) = — "Vg2(Ai, A2), is of the Stackel form and the mechanical system is Hamilton- 
Jacobi separable in the variables Ai, A2. This is the reason of our choice of (3: to cope 
with a new, previously unknown, integrable mechanical system belonging to the class of 
the Neumann problem [25]. 

The solutions of this mechanical system complying with the asymptotic conditions 

Jim ^M. = Hm ^M = , lim (Ai(x), A 2 (x)) = (a 2 , 8 2 ) . (3.2) 

forced by (2.8), will provide all the domain wall solutions with finite tension of the field 
theoretical model. 
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In sphero-conical coordinates the wall tension is: 

fi(Ai - « - srhi / d * [^ (Ai - A *> ) 2 + l^ x - x *{% 2+ ^ <A " A2) 

If a solution, W(Ai, A 2 ), of the PDE: 



is known, the Bogomolnyi arrangement [31] of the wall tension: 

1 f 1 ^ c^A c^Vl^ "\ ^ 1 /* ^ 

n(Ai,A 2 ) = -2 — 2 / dx o yZda ( -r 1 ) + ~2 — 2 / 



<9A ? - 



shows that the absolute minima of the wall tension functional are the solutions of the 
first-order ODE system 

_ =9 n (AiA2) _ i _ =922(AliA2) _ , (3 . 4) 

complying with the asymptotic conditions (3.2). 

However, the PDE (3.3) is no more than the time-independent Hamilton-Jacobi equa- 
tion of the analogous mechanical problem (for zero mechanical energy), such that W(Ai, A2) 
is the Hamilton's characteristic function. Separability implies the existence of solutions of 
the form: 

W(X U X 2 ) = Wi(Ai) + W 2 (X 2 ) 

Integration of the Hamilton-Jacobi equation in general involves hyper-elliptical integrals, 
but the quadratures reduce to simple irrational integrals for the solutions of finite (me- 
chanical) action in infinite (mechanical) time that provide finite tension domain walls. 

The mechanical system is completely integrable in the Arnold-Liouville sense, admit- 
ting two integrals of motion in involution: 

h = \g" A + \g 22 A - 2{x f_ Xl) ((- 2 " Ai)(^ 2 - A1) 2 + (A 2 - a 2 )(A 2 - S 2 ) 2 ) 

h = ^ U A 2 vr 2 + \ g 22 Ax vr 2 - ^ (AiA 2 (A x + A 2 - a 2 - 2S 2 ) - a 2 S 4 ) . 

The finite tension conditions (3.2) require that I\ = I 2 = and the first-order PDE (3.3) 
becomes the ODE system: 



dXi { ' 2 ^1(1 - Ai) ' d\ 2 1 ' 2 ^2(1 -A 2 ) 



(3.5) 



The solution of (3.5) 



W{M, A 2 ) = Li^! ( ^(l-AO + (1 - 25 2 ) arctan ^ 



+ ( iy*R 3 ^ A2(1 _ A2) + (1 _ 2(52)arctan yi^^ ^ ei,e 2 = 0,l 
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is the complete integral of the HJ equation (3.3) such that the ODE system (3.4) reads 

= -R(-ip 2(a 2 -A 1 )(^-A 1 )VA 1 (l-A 1 ) 
dx Ai — A2 

^2 = p j ir 2(a 2 -A 2 )( ( 5 2 -A 2 ) v /A 2 (l-A 2 ) 

A2 — Ai 

encoding all the domain wall solutions as the separatrix trajectories between bounded and 
unbounded motions in the analogous mechanical system. Note that the flow induced by 
the gradient of W is: 

d\i (a 2 - Ai)(<5 2 - Ai)y / Ai(l - Ai) 

i.e., the flow is only undefined at the vacuum (Ai = a 2 ,A2 = 5 2 ), and the focus (Ai = 
a 2 , A2 = a 2 ), points where infinite orbits meet or cross. 

To end this subsection we remark that we must bear in mind two different zones in 
the (Ai, A2)-plane in order to search for the most general wall solutions. The two zones 
are delimited by the straight line A2 = S 2 . In the original variables this line is given 
by the intersection between the conical surface |j + jrz^i — xzf2 — and the sphere 
4>\ + 02 + 03 = i? 2 , which determines two ellipses described on each of the hemispheres. 
We shall refer to these curves as tropical ellipses, using a geographical analogy in the 
sphero-conical coordinates. These so-called tropical ellipses divide the sphere into three 
regions, the intertropical zone (characterized by A2 > S 2 in the sphero-conical plane) and 
the North Polar zone, or "arctic" region, and the South Polar zone, or "ant-arctic" region, 
(characterized by A2 < S 2 in the sphero-conical plane). 

3.3 One more basic wall: Trial orbit in P2 

The orbit A2 = 5 2 , see Figure 8 (left), solves the equation (3.7) and the equation (3.6) on 
this orbit can readily be integrated. We find domain walls that live on the tropical ellipses 



vTW/ \ °" 2 sm h 2 [Rcrcr(x — Xo 

cosh 2 [Raa(x — xq)] — a 2 



Aro»o= v; ; „ -o ^r = s 2 , (3.8) 



which we will refer to as Tropical Domain Walls (TW). It should be understood that 
(3.8) solves (3.6) with e\ = 1 between x — —00 and x = #0, where Ai is decreasing down 
the Ai = axis, whereas it is the solution of (3.6) with e\ = between x = xo and x = 00, 
where Ai increases up the vacuum value a 2 . This kind of domain wall asymptotically 
connects the v 1 and v 2 vacua on the Northern Hemisphere and the v 3 and v 4 vacua on the 
Southern Hemisphere. 

The tension of these domain walls, however, saturates the Bogomolnyi bound 6 : 

ft(TW) = 9 2 9 \W(a 2 ,5 2 ) - W(0,5 2 )\ = 9 R 9 (aa - (1 - 25 2 )arccosa) . (3.9) 
ai — ai 1 1 ai — ai 



6 Although it seems that the tension of these walls depends not only on the value of W at the vacuum 
but also on the value of W(0,5 2 ) at xo, (3.9) is a topological bound. The wall orbit flow is not undefined 
at (Ai = 0, A2 = S 2 ). This point will be clarified further in Section §. 5 . 
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Figure 8: The TW orbit in the P2 rectangle displayed as a solid red line (left) and the TW tension 
density (right) 



The inverse image of (3.8) in § 2 



</>rwi) 



(-l) Kl J R5sinh [Raa{x - x Q )\ 




(-1)^RVS 2 - a 2 



y cosh 2 [Raa(x — xq)] — a 2 
(-l) K3 Mcosh [Raa(x - x )] 




2 [Raa(x — xq)] — a 2 



«1,«2,«3 = 0,1 



(3.10) 



cosh 2 [Raa(x — xq)] — a 2 



helps to elucidate the character of these eight new topological walls, see Figure 9. We stress 
that: (1) The three field components are different from zero. (2) If ^3 = 0, the topological 
wall belongs to C12 (^i = 1) or to C21 (^1 = 0). (3) If ^3 = 1, the topological wall belongs 
to C43 («i = 1) or to C34 (k,i = 0). (4) The sign of ( — l)^ 2 determines the face in § 2 chosen 
by the wall orbit. 



One can check that fi(PMW) > ft(TW) if S > a. Therefore, the stable topological 
walls in these sectors are these latter ones (with 02 7^ 0). This statement is also clear if 
one compares the tension density shown in Figure 8 (right) with the density of the other 
walls living in the same sector, see Figure 10. 

A better understanding of this point comes from the consideration of the wall orbits in 
the meridian 02 = in sphero-conical coordinates. The walls in the tropical zone A2 > 5 2 (m 




Figure 9: Topological wall profiles (3.10) (left) The orbits (right). 
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the C14 and C23 sectors) correspond to the orbit: Ai = a 2 . The analytic solution of (3.7) is: 

S 2 cosh 2 [RSS(x - xq)] 



^TMW 



(x) = a 2 



aTMW 
A 2 



(*) = 



cosh 2 [R88(x - x )] - 8 2 



(3.11) 



Again, (3.11) is the solution of (3.7) if 62 = 
between x = — 00 and x = xq : when A2 increases. 
Between x = xo and x = 00 the choice in (3.7) 
must be 62 = 1 because, then, A2 decreases. 

The orbit starts at the vacuum, reaches the 
equator vertically and bounces back to the vac- 
uum in the sphero-conical plane, see Figure 11 (left). 
The inverse image of this orbit gives the four topo- 
logical walls in C14, C41, C32, C23, Figure 11 (right). 
The wall tension (already given in sub-Section §. 
3.1) for this kind of solution can be computed in 
this context as follows: 




Figure 10: Comparison between the ten- 
sion densities of the two types of walls in 
the same topological sector. 



fi(TMW) 



\W((a 2 ,5 2 ) - W(a 2 ,l) 



R 3 



[c«+(l-25 2 )arccos5] 



(3.12) 

and the same cautionary remarks concerning (3.9) are applicable to (3.12). 

In the polar zone A2 < S 2 we can reproduce the PMW domain wall as a piecewise 
solution in the sphero-conical plane as follows. There is a solution of (3.7) on the orbit 
Ai = a 2 (which is solution of (3.6)): 



^PMWvf _ ^2 ^PMWvf^ 



5 2 sinh 2 [RSS(x-x )] 
cosh 2 [RSS(x — Xq)] — 5 2 



(3.13) 



The orbit also starts from the vacuum and reaches the focus F at the point Xf determined 
by: 



Xf 



1 aS 
xq + arcsmn 



RSS 



VS 2 - a 2 



f V 

I 



ir- 



Figure 11: The orbits of the TMW (2.11) domain walls (left) and the two-stage orbit of the PMW 
(2.10) domain walls (right) displayed in P 2 as a solid red line. 
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Thus, (3.13) is the solution of (3.7), with €2 = 1 between x = —00 and x = —(xf — xo), 

whereas 62 = must be chosen between 3y X j- X Q and x = +00. After hitting the 

focus for the first time, the orbit changes to proceed along the A2 = o 2 axis. The analytic 
solution of (3.6) 

\™»**(X)= ^ 2 [f*>-*0)] , A PMW FP(x) = - 2 (314) 

cosh \R66[x — xo)J — 6 2 

departs from the focus at x = — (xf — xq), solving (3.6) with e\ = 1, rebounds at the pole 
P when x — xo, and travels back with e\ — to reach the focus at x — Xf — xq. There is 
accordingly a continuous sewing of the two stages of the orbit at the focal point, see Figure 
11 (left). 

The inverse map from (3. 13)- (3. 14) to the original field variables produces the four 
topological walls shown in Figure 11 (right) that belong to the sectors C12, C21, C34, and 
C43 of the configuration space. This continuous gluing of two solutions of two different 
first-order equation at a point where the flow is undefined is the bona fide solution of the 
second-order equations, already found in the sub-Section §. 2.2 by direct integration. The 
tension of these walls must be computed in two steps (confirming the result in §2.2): 

f](PMW) = O(PMWvf) + ^(PMW F p) = 

= (\W(a 2 ,6 2 ) - W(d*,d*)\ + \W{a\a 2 ) - W(0,a 2 )|) 

R 3 

— — ~ ?y (66 — (1 — 26 2 ) arccos5) 

4. Composite non-topological and topological domain walls 

4.1 Degenerate families of polar zone non-topological domain walls 

In order to search for more general domain wall solutions we first consider the polar zone: 
A2 E (a 2 ,5 2 ). The first-order equations (3.6)-(3.7) written in differential form: 

(-^^ + _ (-i)^ (4 .1) 



2i?(a 2 - Ai)(<5 2 - AiVAi(l-Ai) 2i?(a 2 - A 2 )(<5 2 - A 2 ) ^2(1 - A 2 ) 

(-l) ei AidAi (-l) e2 A 2 dA 2 

2i?(<7 2 - Ax)^ 2 - Ai)v/Ai(l-Ai) 2i?(a 2 - A 2 )(<5 2 - A 2 )v/A 2 (l - A 2 ) 



-dx(4.2) 



lead to the equation for the mechanical orbit -integrating (4.1)- and the rule for the me- 
chanical time schedule -integrating (4.2)- . Therefore, we obtain a two-parametric (the two 
integration constants Cq, C\ G R) family of domain wall solutions in the implicit form: 



+ f = 2R Cl 



I (a2-Ai)(<y2-Ai)v/Ai(l-Ai) ' J (<t 2 - A 2 )(5 2 - A 2 VA 2 (1 - A 2 ) 

r (-i) £i AirfAi r (-irx 2 dx 2 = 2R(x c) 

J (a 2 -Ai)(<J 2 -Ai)v/Ai(l-Ai) J (a 2 -A 2 )(<5 2 -A 2 )VA 2 (l-A 2 ) " 
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To achieve explicit expressions we instead perform the Euler change of variables: 



1-Ai /l- A 2 „ 2 52 



Si = J " X , 5 2 = y ; < CT2 = ^ < 5 2 < ^1 = ^2 < 5 1 < + °° 

Note the inequalities bounding the new variables Si,i = 1,2, in terms of the old param- 
eters a and 5 as separation constants. A cosmetic change of notation in favor of <J\,<J2 
is also introduced to make the formulas more symmetric. The quadratures (4.1)-(4.2) are 
rationalized: 

1 y f (-1)^(1 + s 2 )ds t 1 A f (-lY'dSj _ 

6 2 a 2 R J (a 2 - s 2 )(a 2 - s 2 ) ° X ' 8 2 a 2 R (a 2 - s 2 )(a 2 - sj) °° " 

(4.3) 

Plugging the simple fraction decompositions 



2 



1 S 2 a 



2 



("? - '?) K 2 - «?) ^ 2 - * 2 V of - o\ - s 2 
1 + s 2 1 / S 2 a 2 



(*l-s 2 )(4-s 2 ) a 2 -6 2 \a 2 ~s 2 a 2 -s 2 
in (4.3) further simplifies the quadratures: 

E / =S 2 R(x-C -a 2 C 1 ), J] f i - Z i^ = ° 2 R(x-Co-6 2 C 1 ) . 

i=l J °1 S i i=l J ^ S i 

(4.4) 

Integration of the ODE's (4.4) provides the two-parametric family of domain wall solutions: 

arccoth — + arccoth — = R 8 6 (x — xq) , arccoth — + arctanh — = R a a (x — xq + C) > 
(72 (72 <r\ a 1 

(4.5) 

now in terms of other two integration constants: xo = Co + cr 2 Ci, ( = (a 2 — 5 2 )Ci, 
xo, C £ ^- Setting the value of the integration constant xq fixes the "center of mass" of the 
wall, whereas different values of £ determine the different wall orbits uniquely. We have 
re-defined the Si variables in the form si —> ( — l) ei si, S2 — > (— l) e2 <§2, taking advantage of 
the parity properties of the inverse hyperbolic functions. 
The addition formulas for the hyperbolic functions: 

tanh (arccoth + arccoth q) 



1 + pq tanh (arccoth p + arctanh q) 
allow us to invert (4.5) and we find: 
a 2 {si + s 2 ) tanh r R § - _ , ^ h ox(o\ + 8182) = tanh [R(J - {x _ Xq + c)] _ t2 

0-2(^2 + S 1 S 2) Oi(si + S 2 ) 

(4.6) 

This is an algebraic linear system in the "Vieta variables" A = si + 52 , B — s\S2'. 

a 2 A - t\B t\o\ , ait 2 A - B = a\ , 
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solvable by means of Cramer's rule: 

(S 2 -a 2 ) tl 



A(x;x X) 



a5 (aS — 6atit2) 



si, 52 are by definition the roots of the quadratic equation s 2 — As + B — 0: 



, N A(x) + JA 2 (x) -4B(x) , N A(x) - JA(x) 2 - AB(x) 
s i( x ) = g ' S2 ^ = 2 ' 

and thus we find explicit expressions for the family of domain wall solutions of (3.6) and 
(3.7) in the polar zones: 

aPZW- 



Ar ZW (*;* ,C) 



1 + sf(x;x ,Q 



1 + sf,{x;x ,Q 



(4.7) 









/ Mx \ 




-4 / -1 

-4? 


2 \ 4 



Figure 12: Graphics of the domain wall components (4.8, 4.9, 4.10) for: (left) Xo = C = 0> an d 
(right) x = 1, C = 3. 

although we stress that in the formula (4.7) the signs of (3.6) and (3.7) must be chosen to 
fit with the stages where Ai and A2 are respectively increasing or decreasing. 

We return to Cartesian coordinates in R 3 using formula (3.1). (f> 2 , (p 2 and <fi 2 are given 
in terms of s\ + s 2 and s 2 s 2 . The analytic expressions depending on the spatial coordinate 
x, and the integration constants 71 and 72, are: 

(-1) K1 R5 (Sa-Sahh) 



tf zw (x;x ,C) 



^ ZW (*;*o,C) 



^ zw (*;*o,C) 



(<S 2 - a 2 ) 2 1\ + a 2 a 2 t 2 t 2 - 255aat 1 t 2 + 5 2 S 2 

(-i)*2R(p-<T*)t ly /r-% 

(S 2 - a 2 ) 2 1\ + o 2 aH\t 2 - 2SSaaht 2 + S 2 P 
(-l) K3 R6(Sa - Saht 2 ) 



(4.8) 
(4.9) 
(4.10) 



'(S 2 - a 2 ) 2 1 2 + o 2 oH\t\ - 2SSaaht 2 + S 2 S 2 
The meaning of the integration constants xq, £ is now clear: 

1. Besides fixing the center of mass of the composite wall, xq sets the point xj = xq, where 
the field profiles touch the foci. At this point t\ = and : 

(^ zw ) 2 (*o;* ,C) - R2 ° 2 > (^ ZW ) 2 (^o;xo,C) = , (^ zw ) 2 ( W o,C) - i?V . 

2. £ determines the orbit of the non-topological wall and fixes the relative coordinate 
between the two centers of these composite walls. 
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In Figure 12 the graphics of the three components 
of the non-topological domain wall field profiles for two 
sets of values of (xo,0 are snown: ^i ZW (blue lines) 
tend to -5, $ zw (red lines) tend to 0, and 
(brown lines) tend to 5 at x — » +oc and x — oc; 
i.e., the wall solutions go to the same vacuum at both 
ends of the straight line, which determines the non- 
topological character of these solutions that belong to 
the Cii, C22, C33, and C44 sectors. In the plots of the 
wall orbits in § 2 drawn in Figure 13, it is seen that 
all the orbits start and end at the same vacuum and 
cross through the opposite focus. More interestingly, 
the tension densities of these walls depicted in Figure 
14 unveil their character as composite extended ob- 
jects: they are non-linear superpositions of the PMW 
and TW basic walls. 




Figure 13: Orbits in S 2 for three 
PZW domain wall solutions: (1) 
(71 = 0,72 = 0), red. (2) ( 7l = 
-1,72 = 1), brown. (3) (71 
1,72 = 6), blue, (right). 




Figure 14: Non-topological wall tension densities for:(l) (xq = 0, £ = 0), red. (2) (xq = 1,£ = 2), 
brown. (3) (xq = — 1,C = 5), blue. 



To further explain these statements, we discuss the orbits in the P2 rectangle. In 
Figure 15 (left), the plots of three of these wall orbits are shown. The (xq = 0, ( = 0) orbit 
(green line) starts at the vacuum, hits the Ai = 0-axis at the mid-point between A2 = S 2 
and A2 = a 2 , and runs to the focus F. Then, the orbit returns back to the vacuum through 
the same path in reversed sense. 

The (xo = 1, C — 2) orbit (brown line) starts at the vacuum point, hits the Ai = 0-axis 
at a point below A2 = ~^ , and run to the focus. Then, the orbit runs again to intersect 
the Ai = 0-axis at the symmetric point over A2 = ~^ ? and travels back to reach the 
vacuum. 

The (xq = — 1,C — 5) orbit (black line) follows a similar pattern in an extreme way. 
The trajectory starts from the vacuum, running almost parallel to the A2 = S 2 orbit, hits 
the Ai = 0-axis, and returns slightly departing from the same way until, close again to the 
vacuum, it turns down toward the focus, almost parallel to the Ai = a 2 -axis. Just at the 
focus, the orbit turns to the left, again hitting the Ai = 0-axis at a point extremely close 
to the North Pole where it turns back almost parallel to the Ai = a 2 -axis. 

Again approaching the focus, the orbit finally turns up almost vertically to end at the 
vacuum. Needless to say, the follow-up of the orbits described here can be interpreted in 
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the inverse sense due to the symmetry x — >> — x of the system. A very important point to 
emphasize is that the focus F is a conjugate point with respect to the vacuum V: a point 
where all the orbits of a congruence starting in V pass through. This will have consequences 
in the stability of these domain walls, an issue to be analyzed in a forthcoming publication. 

The (xq = — 1,C = 5)-orbit is particularly interesting: it is very close to the gluing 
of the two trial orbits in the C\2 and C34 sectors of the system. This confirms that these 
new domain wall solutions are nonlinear superposition of two basic walls. The Bogomolnyi 
trick provides the tension of any wall in this family : 



ft(PZW) = 



at 



at 



+- 



at 



at 



|W(a 2 ,5 2 )- W(0,5 2 )| + 

(\W(a 2 ,a 2 ) - W(0,a 2 )\ + \W(a 2 ,5 2 ) - W(a 2 ,a 2 )\) 



R 3 



at 



an 



7T 



a a + SS + (1 — 26 )(arcsina — — — arcsin5) 



(4.11) 



= fi(PMW) + fi(TW) 



The formula (4.11) is a remarkably result. First, it means that the tension of all the walls 
in the family is the same. Second, the tension is equal to the sum of the tensions of the two 
basic walls that live in the same topological sectors. This wall tension sum rule is another 
confirmation that the polar domain walls are composed of two basic domain walls. 





Figure 15: Graphics of the domain wall orbits in P2 displayed as solid lines: Polar zone non- 
topological domain walls (left) and Tropical zone topological walls (right). 



4.2 Degenerate families of tropical zone topological domain walls 

If 5 2 < A2 < 1, the tropical zone, the inequalities satisfied by the Si variables are: 

al <a\ < s\< +oc , < s\ < g\ 

Integration of the ODE's (4.4) gives the domain wall solutions in this zone 

arccoth — + arctanh — = R S 5 (x — xq) , arccoth — + arctanh — = R a a (x — xq + () 
&2 &2 &i o\ 



(4.12) 
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The only difference with respect to (4.5) is that, s\ also being smaller than cr^, there are 
two arcth functions entering the solution. The subsequent linear system in Vieta variables 
becomes 

a 2 t 2 A - B = &% j (J\t\A -B = a\ . 



Cramer's rule dictates the solutions: 



A(x;xqX) 



(S 2 



(7 



S (cra5t\ — 8a 2 t2) 



a (aSSt2 — crS 2 ti) 
5 (aa5t\ — 5a 2 t2) 



si(x) 



A(x) + y/A 2 (x) -4B(x) 



s 2 (x) 



A(x) - ^A(x) 2 - 4B(x) 



and we have a new two-parametric family of tropical domain walls of equations (3.6) and 
(3.7) with the appropriate signs: 



AP W (*;*o,C) 



1 



, Aj zw (x;x ,C) 



1 + sf(x;x ,Q 



1 



1 + sl{x;x Q ,Q 



(4.13) 





Figure 16: Graphics of the tropical zone topological wall components (4.14, 4.15, 4.16) for: (left) 
.To = 0, C = 0, and (right) xq = 1, ( = 3. 

Back in Cartesian coordinates in field space, we find the new family of domain wall 
solutions: 



APW(x;x„,C) = jg^j (4,4) 
yj (vvt x - SSt 2 ) + (S 2 - a 2 ) 2 

a tzw (x;xo>c) = (-irw-f)tt (4 , 5) 

yj (aah - 88t 2 ) 2 + (<5 2 - a 2 ) 2 
y/ (vat! - 55~t 2 f + (5 2 - a 2 ) 2 



(4.16) 



Figure 16 plots the field profiles of two domain walls in this family. They run asymptotically 
into different vacua and are thus topological solutions living in topological sectors of the 
configuration space. In general, the topological wall solutions (4.14, 4.15, 4.16) belong to 
the topological sectors: C24/C42 and C13/C31. Their wall tension densities have the form 
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Figure 17: Topological wall tension densities for:(l) (xo = 0, ( = 0), red. (2) (xq 
brown. (3) (xq = — 1,C = — 7), blue. 



-1,C = -2), 



shown in Figure 17, again suggesting the composition of two basic walls: TMW and TW 
in this case. 

The integration constants xo and £ determine the center of mass and relative coordi- 
nates of the composite walls, just as in the non-topological polar walls. 

The orbits in S 2 of three topological walls 
are plotted in Figure 18. They join anti-podal 
vacua and there are no conjugate points in these 
congruences, a fact that offers a strong hint 
of stability from the Morse index theorem, see 
[32]. 

Finally, we describe the topological wall or- 
bits in the P2-rectangle in Figure 15 (right). 
Unlike the non-topological wall orbits that are 
formed by six, their orbits are composed of four 
stages. The change of stage takes place when 
the orbit hits either the A2 = 1 (the equator) 
or the Ai = (the <\>\ = meridian) edges. 
The topological wall orbits do not pass through 
the foci, the points where the flow is undefined, 
and therefore their tension is a true Bogomolny 
bound, depending only on the values of the fields at the vacua: 




Figure 18: Orbits in S 2 of three topological 
domain wall solutions: (1) (xo = 0, £ = 0), 
red. (2) (x = -1,C = -2), brown. (3) 
(a? = -l,C = -7), blue. 



Q (TZW) = 



R 3 ' 

2 2 

af — «3 



(\W(a 2 ,5 2 ) - W(0,S 2 )\ + \W(a 2 ,S 2 ) - W(a 2 ,l)\) 
{a a + 55 + (1 — 25 2 )(arcsina — arcsin^)) 



(4.17) 



= fi(TMW) + ft(TW) . 



Accordingly, the topological walls (4.14, 4.15, 4.16) form a degenerate family of BPS walls 
in tension, confirming the stability of these topological defects. 

5. Further comments 

In [30], we discussed the stability properties and quantum features of the domain walls 
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discovered in [15]. We plan a similar analysis of the domain walls described in this paper 
in future research. Nevertheless, we shall briefly comment on these points in this last 
Section: 

1. Stability. It is compelling to wonder about the stability of these domain walls. The 
answer is as follows: 

• The topological basic walls TMW as well as the topological TW (and their anti's) 
are stable . 

• The topological basic walls PMW (and their anti's) are unstable . 

• The composite non-topological walls are unstable . 

• The composite topological walls are stable . 

The arguments to support these claims are based on: 1) The saturation of Bogomolny 
bounds, see [31]. 2) The application of the Morse index theorem, see [32]- [33]. 3) The 
computation of Jacobi fields, ( not given in this paper), see [34]. 

2. One-loop wall tension shift. Following the work in [35] on the supersymmetric kink we 
computed the one-loop mass shift to the only stable kink found in [15] in Reference [36]. 
This work has been extended to the massive model with target space § 3 in [37]. In this case 
there are two stable topological kinks but the strategy used was the spectral zeta function 
regularization developed in the papers [38]- [39]- [40] on linear sigma models with several 
scalar fields. 

Finally, a quick remark about what happen if the parameter 6 took other values. 
First, if 6 > 1 the four ground states become imaginary but the maxima at the intersection 
between the Equator and the 02 = 0-meridian become minima. There would be only two 
vacua on these two points and the structure of the moduli space of domain walls would 
be very similar to the structure described in [15], although the analytic expressions of the 
domain wall solutions would differ. Second, if 5 = a the ground state and the foci would 
coincide and the two types of wall orbits -connecting either antipodal or non-antipodal 
pairs of vacuum points- would be topological. All the domain defects would be stable 
composite topological walls in this case, whereas the basic topological walls would all be 
stable. 
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